Ramsey interferometry with atoms and molecules: two-body versus many-body phenomena 
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We discuss the frequency and visibihty of atom-molecule Ramsey fringes observed in recent experiments 
by Claussen et al. [Phys. Rev. A 67, 060701 (2003)]. In these experiments a *^Rb Bose-Einstein condensate 
was exposed to a sequence of magnetic field pulses on the high field side of the 155 G Feshbach resonance. 
The observed oscillation frequencies largely agree with the theoretically predicted magnetic field dependence 
of the binding energy of the highest excited diatomic vibrational state, except for a small region very close to 
the singularity of the scattering length. Our analytic treatment of the experiment, as well as our dynamical 
simulations, follow the magnitude of the measured oscillation frequencies as well as the visibilities of the Ram- 
sey fringes. We show that significant deviations from a purely binary dynamics, with an associated binding 
frequency, occur when the spatial extent of the molecular wave function becomes comparable with the mean 
distance between the atoms in the dilute gas. The experiments thus clearly identify the conditions under which 
diatomic molecules may be identified as a separate entity of the gas or, conversely, when the concept of binary 
physics in a many-body environment is bound to break down. 

PACS numbers: 03.75.Kk, 34.50.-s, 05.30.-d 



I. INTRODUCTION 

The study of ultra-cold molecules produced in atomic Bose 
and Fermi gases is one of the most important developments in 
cold atom physics Qj. One of the most successful experimen- 
tal techniques to produce molecules in atomic Bose-Einstein 
condensates utilised adiabatic sweeps of a magnetic field tun- 
able Feshbach resonance level across the zero energy thresh- 
old of the colliding atoms, during which a highly excited di- 
atomic vibrational bound state could be efficiently populated 
QmQEll^. In this paper we study a coherent superposition 
of an atomic condensate and molecules. The production of 
such a coherent superposition in a degenerate gas of bosons, 
demonstrated in Refs. |7, is a remarkable achievement, 
given the marked difference in the relevant energy and time 
scales for condensates and molecules. These experiments em- 
ployed a sequence of two magnetic field pulses each of which 
rapidly approached the position of the zero energy resonance 
(the magnetic field strength at which the scattering length has 
a singularity) on the side supporting the highly excited di- 
atomic vibrational bound state. The two pulses were sepa- 
rated by an evolution period of variable duration as a function 
of which the oscillations in the final condensate populations 
were observed. This achievement offers the possibility of pre- 
cise interferometric measurements of the energies of the rele- 
vant states through the measurement of bulk properties of the 
gas. It is the theory of such measurements we present in this 
paper. 

We set up the theory in a way that reflects properly the 
analogy with Ramsey - separated pulse - interferometers with 
the splitting and recombination occurring during the magnetic 
field pulses at the beginning and the end of the magnetic field 
sequence, respectively. We give a complementary account of 
the experiment which emphasises the interferometer picture 
and will, we believe, be of greater utility to those analysing 
experiments than just the simulation on its own can do. 

To provide a simple intuitive physical description of the ex- 



periments, we first consider the evolution of a single pair of 
cold trapped atoms subjected to the experimental pulse se- 
quence. We show that the first pulse splits the initial state of 
the pair into three different components (interferometer arms). 
The bound state component evolves most quickly. The other 
components are a superposition of excited quasi continuum 
states and the lowest energetic quasi continuum state, which 
are the two-body analogues of the experimentally observed 
13] burst of atoms and the remnant Bose-Einstein condensate, 
respectively. After a variable period of time, referred to as 
the evolution period, during which the bound and free states 
evolve almost independently, the second pulse recombines the 
components. The final state of the pair depends sensitively on 
the phase differences built up between the amplitudes of the 
various components during the evolution period. As a con- 
sequence, oscillations in the resulting populations of the three 
components of the pair wave function are observed as the time 
between the pulses is varied. We shall show that the oscil- 
lation frequency is determined by the energy of the molec- 
ular bound state populated in the process. The visibility of 
these oscillations, i.e. Ramsey fringes, depends on both the 
efficiency of populating the bound state during the first pulse 
and of dissociating the molecules during the second pulse of 
the sequence. We observe that these quantities are very sen- 
sitive to both the details of the underlying two-body physics 
and the precise duration and shape of the pulses. We provide 
analytic estimates of the molecular production and dissocia- 
tion efficiencies and give an analogy between the function of 
the pulses and the concept of two colour photo-association of 
molecules. 

We extend this description of a single pair of atoms in or- 
der to apply it to a gas of atoms by using a coarse graining 
procedure. In this approach the populations of the different 
components of the gas are obtained, following classical prob- 
ability theory, by averaging over the microscopic transition 
probabilities for all pairs. We show that this straightforward 
but systematic extension of the two-body physics recovers the 
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experimental observations rather well, in terms of the oscilla- 
tion frequency, the amplitude and the relative populations of 
all three components present in the gas. 

In order to disentangle the two-body and many-body as- 
pects of Ramsey interferometry with cold gases, we have also 
used the genuinely many-body microscopic quantum dynam- 
ics approach to the dynamics of partly condensed Bose gases 
II9I [Toll . We base our full simulations on the non-Markovian 
non-linear Schrodinger equation of the first order microscopic 
quantum dynamics approach. This equation describes the evo- 
lution of the condensate along with the dynamics of pair corre- 
lations. The method describes fully and quantitatively the evo- 
lution of the states produced by magnetic pulse shape induced 
rapid changes in the Hamiltonian describing the inter-atomic 
potential. In particular, the microscopic quantum dynamics 
approach incorporates the entire binary collision dynamics in 
a non-perturbative manner 

We shall see that the many -body phenomena affect both the 
frequency and visibility of the Ramsey fringes. In the first 
step in this direction, we have found an analytical solution 
of the non-Markovian equation using a perturbative approxi- 
mation. In agreement with the experimental observations, we 
show that the frequencies of the Ramsey fringes are shifted 
upwards with respect to the value associated with the energy 
of the weakly bound molecular state populated in the process. 
This shift persists for all values of the magnetic field, it does, 
however, become hard to discern above Bevoive (the value of 
the magnetic field strength between the pulses) higher than 
roughly 158 G. Away from the zero energy resonance the shift 
is given, as one might reasonably expect, by roughly twice the 
chemical potential. The shift is typically a small correction to 
the oscillation frequency, and a time scale much longer than 
the typical pulse durations of the experiments |01 would be 
required for it to be fully resolved. The visibility of the Ram- 
sey fringes is reduced in comparison with the two-body result 
due to the non-linear dynamics of the atomic condensate mean 
field. 

Our analytic studies apply when the experiments 0, @1 are 
operated in what we shall term the "interferometer mode", 
i.e. when the magnetic field strength Bevoive during the evo- 
lution period is sufficiently far away from the zero energy res- 
onance for two- and many-body phenomena to be disentan- 
gled. In this regime of magnetic field strengths Bevoive the gas 
is weakly interacting during the evolution period and no sig- 
nificant population exchange between the bound and unbound 
components occurs. The different components are then or- 
thogonal as they should be in an ideal Ramsey interferometer. 
Claussen et al. jH also performed measurements in the close 
vicinity of the zero energy resonance, where the gas is subject 
to a non-equilibrium dynamics even between the pulses. For 
all experimentally studied cases we have performed full sim- 
ulations of the non-Markovian equation for the trap geometry 
and atom number employed in the experiments. In the regime 
very close to the zero energy resonance the non-linear aspects 
of the dynamics are especially significant, and the description 
of the experiment in terms of an ideal Ramsey interferometer 
is not as clear, because the three components of the gas are 
no longer orthogonal. The full numerical simulation does not 



rely on this assumption and is, therefore, essential in the com- 
parison with experimental data for this regime. We should 
emphasise that in this region the idea that we can clearly sep- 
arate molecules and atoms is also breaking down as the spatial 
extent of the molecular bound state wave function approaches 
the mean distance between the atoms in the gas j ll]. This 
makes the simple picture invalid as well as any theory that 
relies on such a distinction. 

One might be misled in this analysis if one focused purely 
on the component of the molecular bound state in the bare 
resonance level that remains, of course, very small in extent. 
Our analysis shows, for instance, that the frequently used two 
level mean field approach of Refs. (ll[ll[ll[ll[!l[ll[ll, 
fails entirely to describe the interferometric experiments of 
Refs. 0,01, despite the fact that it can provide accurate predic- 
tions on the final molecular fractions in Feshbach resonance 
crossmg experiments In Ref. j2lll a modified two 

level mean field approach was applied to the observations of 
Ref. m, in which the energy of the bare Feshbach resonance 
level was tuned to recover the binding energy of the molecules 
produced in the experiments Q |^ . This approach recovered 
the magnitude of the fringe frequencies of Ref. 01 after ad- 
justing the densities used in the theory. We believe, however, 
that such a treatment is inconsistent because any two level ap- 
proach crucially lacks the excited energy levels of the back- 
ground scattering continuum I20I1 . This treatment obviously 
leads to the absence of the experimentally observed burst com- 
ponent of atoms 131 in the theory. The two level approach of 
Ref. ^IJ also inevitably misrepresents the molecular bound 
state I22II produced in the experiments jH, as this bound state 
is even dominated by its component in the background scat- 
tering continuum. 

The Ramsey fringes predicted by our dynamical simula- 
tions of the non-Markovian equation recover all the system- 
atic trends observed in the experiment. For magnetic fields 
away from the zero energy resonance, the Ramsey fringes are 
essentially undamped on the comparatively short time scales 
of the pulse sequences and their frequencies are only slightly 
shifted with respect to the value determined from the energy of 
the weakly bound molecular state. Closer to the zero energy 
resonance, a significant decoherence of the Ramsey fringes 
has been observed in the experiment, a feature also predicted 
by our many-body approach. We also recover the substantial 
upward shift of the fringe frequencies for the magnetic field 
strengths closest to the zero energy resonance, for which the 
shift was attributed in Ref. I^ Esil . We show that this shift can 
be interpreted in terms of a measure of the extent to which 
diatomic molecules in the highly excited vibrational state can 
be considered as a separate entity of the gas. Our studies also 
strongly indicate that the significant experimentally 0] ob- 
served loss of condensate atoms, which is particularly clear 
for longer pulse sequences with the magnetic field strength 
^evolve in the close vicinity of the zero energy resonance, is 
essentially determined by the same parameters of the binary 
interactions that determine the frequency and visibility of the 
fringe pattern. This observation implies that deeply inelastic 
loss phenomena such as resonance enhanced spin relaxation 
of **^Rb atoms I24I1 . are unlikely to play a significant role in 



3 



the experiments | 

Both the two-body and many-body pictures presented in 
this paper reveal the remarkable sensitivity of these experi- 
ments to the details of the underlying two-body physics, de- 
tails that are clouded over in experiments involving molec- 
ular production using linear crossing of a Feshbach reso- 
nance fl9ll20ll . The interferometric technique, therefore, of- 
fers a rather unique probe of the microscopic collision physics 
through observing the macroscopic bulk properties of the gas. 



II. PHYSICAL ORIGIN OF THE RAMSEY FRINGES 

In this section we provide a physical explanation for the ex- 
perimental fringe patterns reported in Refs. 001 in terms of 
general properties of a Ramsey interferometer. We then pro- 
vide a very direct determination of the frequency and visibil- 
ity of the Ramsey fringes in terms of the dynamics of a single 
pair of atoms. We then present the first order many-body cor- 
rections to the frequency and visibiUty obtained from the mi- 
croscopic quantum dynamics approach of Refs. f^,^. This 
allows us to distinguish pure two-body contributions to the 
interference between atoms and molecules from many-body 
phenomena related to the coherent nature of the partially con- 
densed Bose gas. 



A. Atom-molecule coherence in a Bose-Einstein condensate 

In the experiments of Refs. 0,01 a Bose-Einstein conden- 
sate of '^^Rb atoms in the {F = 2, ruf - -2) hyperfine state 
was exposed to a sequence of spatially homogeneous mag- 
netic field pulses on the high field side of the 155 G Fesh- 
bach resonance. A typical pulse sequence l25ll is depicted in 
Fig. n The pulses, at the beginning (#1) and the end (#2) 
of the sequence, each consisted of linear downward and then 
upward ramps, which were separated by a hold period of con- 
stant magnetic field strength, close to the position of the zero 
energy resonance at Bq = 155.041 G liSJ. The evolution period 
of constant magnetic field strength Bevoive separated the pulses 
by an amount of time fevoive- In the course of the experiments 
0|8||, the final densities of atoms were recorded after expand- 
ing the cloud, and the measurements were repeated with vari- 
able evolution times fevoive and magnetic field strengths Bevoive- 



1. Interference between different components of a partially 
condensed Bose gas 

Donley et al. 0] have identified two different components 
of the final atomic cloud: a remnant Bose-Einstein conden- 
sate and a burst of atoms with a comparatively high mean ki- 
netic energy of roughly Ekin/^B = 150 nK (k^ - 1.3806505 x 
lO^^^* J/K is the Boltzmann constant). The magnitudes of the 
burst and remnant condensate fractions exhibited an oscilla- 
tory dependence on the evolution time fevoive of the magnetic 
field pulse sequence. The existence of a third component of 
missing atoms was inferred from the difference between the 
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FIG. 1: A typical sequence of magnetic field pulses (solid lines) ap- 
plied in the interferometric experiments of Refs. QSl- The magnetic 
field strength at the minimum of the pulses #1 and #2 of B^i„ = 155.5 
G is closest to the zero energy resonance at Bq = 155.041 G 
01 (dashed, horizontal line), and leads to large positive scattering 
lengths on the order of 10000 OBohr (fsohr = 0.052918 nm is the Bohr 
radius). The variable magnetic field strength and duration of the evo- 
lution period of the pulse sequence are denoted by Bevoive and fevoive i 
respectively. 



initial and total final numbers of detectable atoms. Donley et 
al. 0] suggested that the missing atoms indicated molecular 
production, and the oscillations were interpreted in terms of 
an interference between atoms and molecules during the evo- 
lution period of the pulse sequence. 

This interpretation has been confirmed by subsequent the- 
oretical studies fiol |25 EtIi . Figure |2] shows a theoretical 
prediction of the final densities of atoms versus fevoive for the 
pulse sequences and the initial conditions corresponding to 
Fig. 6 of Ref. 0. The predictions in Fig. |2] were obtained 
from the microscopic quantum dynamics approach to dilute 
Bose-Einstein condensates Lid. 28.] . This gives a fringe pat- 
tern which compares rather favourably in both its qualita- 
tive and quantitative features with the observations reported 
in Ref. |7]. The same theoretical approach also recovers the 
experimentally observed kinetic energy of the burst atoms on 
the order of fcRXl50 nK fioll. 



2. The different regimes of molecular vibrational frequencies 

The frequencies of the fringes have been identified as the 
vibrational frequencies of the highest excited diatomic vibra- 
tional bound state 0b(fievoive) at the magnetic field strength 
^evolve of the cvolution period 0, through their remarkable 
agreement with the result of coupled channels calculations 
(cf. Fig. 13. The marked magnetic field dependence of the 
binding energy £'b(Bevoive) associated with the bound state 
0b(fievoive) IS depicted in Fig. |3l We note that the range of 
frequencies measured via bulk properties of the gas extends 
over three orders of magnitude from several kHz to the MHz 
regime. The associated oscillation periods are very much 
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FIG. 2: Predicted numbers of condensate and burst atoms at the end 
of the magnetic field pulse sequence as a function of /evolve for a mag- 
netic field strength of i?cvoivc = 159.84 G in the evolution period of 
the pulse sequence. The total number of atoms is 17100 (horizon- 
tal line). The predictions were obtained from numerical solutions of 
the non Markovian non-linear Schrodinger equation <28> . including 
a fully quantum mechanical description of the inhomogeneity of the 
gas 1231 and a two-channel description of the low energy binary col- 
lision dynamics. The magnetic field pulse sequence (cf. Fig.0, trap 
frequencies, and initial ground state condensate densities correspond 
to the measurements in Fig. 6 of Ref. 01 . 



shorter than typical time scales set by the bulk motion of dilute 
Bose-Einstein condensates, which are typically on the order of 
the trap period of about 1/10 of a second in Refs. 

For the regime of the high MHz oscillation frequencies the 
spatial extent of the molecular wave function is typically on 
the order of several nanometres iHHi . while the length scale 
of the atomic cloud is set by the trap length of about 3000 
nm. In view of these enormous ditferences in the orders of 
magnitude the visibility of the Ramsey fringes is a remark- 
able phenomenon in its own right. In the opposite regime of 
low kHz oscillation frequencies, i.e. very close to the zero en- 
ergy resonance at Bq - 155.041 G, Claussen et al. jSji] have 
observed pronounced upward shifts of the fringe frequencies 
with respect to the vibrational frequencies of the molecules. 
The frequency shifts were accompanied by a strong damping 
of the fringe pattern. This damping suggests an enhancement 
of the influence of many-body phenomena related to the co- 
herent nature of the dilute Bose-Einstein condensate on the 
fringe pattern, as the magnetic field strength approaches the 
zero energy resonance. 



Ramsey interferometry with atoms and molecules in a 
Bose-Einstein condensate 



We shall now identify the different elements of a Ramsey 
interferometer from the pulse sequence in Fig. [2 Our analysis 
will show that in the regime of comparatively high molecular 
vibrational frequencies the atom-molecule oscillation exper- 
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PTG. 3: Binding energy Eb(B) of the highest exited vibrational 
molecular bound state if>b(B) of two ^^Rb atoms in dependence on the 
magnetic field strength B. The figure shows comparisons between 
full coupled channels calculations IM l2^ (circles), the two channel 
approach of Refs. fT^l23l (solid curve), and single channel analytic 
estimates ficll of the binding energy (dotted and dashed curves) in 
the vicinity of the zero energy resonance at Bq = 155.041 G. While 
the dotted curve refers to the universal estimate, just depending on 
the scattering length a, the dashed curve also accounts for the van 
der Waals interaction -C(,/r^, at large inter-atomic distances r, in 
terms of the mean scattering length a fiollscll . The dotted-dashed 
line indicates the linear dependence of the energy £,.es(S), associ- 
ated with the Feshbach resonance level 0res(''). on the magnetic field 
strength. Here and throughout this paper, we have chosen the zero 
of the energy, at each magnetic field strength B, as the energy of two 
asymptotically free *''Rb atoms in the (F = 2, mf = -2) state. 



iment |7|, |8|] can be directly interpreted in terms of Ramsey 
interferometry, while in the opposite low frequency regime a 
loss of coherence is expected to wipe out the fringes. 



1. Interpretation of the magnetic field pulses in terms of beam 
splitters 

We shall first show that the pulses #1 and #2 of Fig.^pro- 
vide the analogue of beam splitters through rapid tuning of 
the inter-atomic forces in the gas between the weakly and the 
strongly interacting regime. This is achieved by the variation 
of the scattering length a{B), which characterises the static bi- 
nary interactions in dilute Bose-Einstein condensates and de- 
pends on the magnetic field strength B in accordance with the 
formula: 



a{B) = flbg 1 



AB 
B-Bo 



(1) 



Here abg - -443 cfBohr (flBohr = 0.052918 nm is the Bohr ra- 
dius) and AB = 10.71 G |8| are the background scattering 
length and the resonance width, respectively. At the min- 
imum of the magnetic field pulses the stationary magnetic 
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field strength of B^in - 155.5 G thus implies a large posi- 
tive scattering length of about 10000 OBohr- which is close to 
the mean inter-atomic distance of about 12000flBohr in the di- 
lute gas. The singularity of the scattering length in Eq. ([Q is 
a decisive feature of all zero energy resonances in ultra-cold 
collision physics and indicates the emergence of a highly ex- 
cited diatomic vibrational bound state (pb{B) at the magnetic 
field strength Bq of the zero energy resonance. This molecular 
bound state exists just on the side of positive scattering lengths 
(cf. Fig. 13, and its mean internuclear distance is largely deter- 
mined by <r> = a{B)/2 

The favourable overlap between the diatomic energy states 
and the many-body state of the gas provides the necessary 
conditions for producing the desired strongly correlated pairs 
of atoms in the molecular bound state 0b(fievoive)- during the 
upward ramp at the end of the first magnetic field pulse. The 
association of bound molecules at the final magnetic field 
strength Bevoive of Pulse #1 is accompanied by the production 
of correlated pairs of atoms in the continuum of positive en- 
ergies above the two-body dissociation threshold. The occu- 
pation of these continuum states provides a first burst compo- 
nent that diffuses during the evolution period. The first mag- 
netic field pulse thus splits the virtually uncorrected initial 
Bose-Einstein condensate into three components; the remnant 
Bose-Einstein condensate, bound molecules and burst atoms. 



2. The evolution period 

In the ideal limit of coherent Ramsey interferometry, 
the gas is weakly interacting during the evolution period, 
i.e. Bevoive is sufficiently far from the zero energy resonance 
that the scattering length a(B) is much smaller than the mean 
distance between the atoms in the dilute gas. The remnant 
condensate, burst and bound molecular components are then 
virtually orthogonal and evolve coherently. The amplitude of 
the molecular component in the state <pb{Bewo\ve) thus acquires 
a phase shift of Aip - |£'b(Bevoive)|revoive/^ (cf. Fig. 13 relative 
to the slowly varying amplitudes of the burst and condensate 
components. 



3. The second magnetic field pulse 

The second magnetic field pulse gets all three components 
to overlap and thereby probes the phase differences between 
their amplitudes. The final densities of the three components 
are therefore sensitive to the phases of their amplitudes at the 
end of the evolution period, which explains the emergence of 
a fringe pattern as a function of fevoive as well as its angular 
oscillation frequency of We - liibCfievoive)!/^- 

We note that in order to observe coherent Ramsey fringes, 
the different components of the system need to be orthogonal 
during the evolution period. In Ramsey interferometers with 
light or atomic beams the beam splitters provide the spatially 
separated, i.e. orthogonal, components. In the atom-molecule 
coherence experiment orthogonality is provided by the ide- 
ally large difference in the orders of magnitude between the 



mean inter-atomic distance in the atomic burst and condensate 
components on one hand and the comparatively small bond 
length of the bound molecules during the evolution period on 
the other hand. In both cases the lack of overlap prevents any 
significant exchange between the components during the evo- 
lution period. 



C. Frequency and visibility of the interference fringes 

In this subsection we determine the visibility and oscillation 
frequency of the Ramsey fringes in the remnant Bose-Einstein 
condensate and provide quantitative corroboration of the more 
physical arguments provided in Subsection lllBI We shall pro- 
vide the main physical ideas of the derivations in terms of the 
dynamics of a single pair of atoms and present the results of 
a genuine many-body treatment to clearly separate two-body 
phenomena from those specifically related to the coherent na- 
ture of the partially Bose condensed gas. The detailed deriva- 
tions on the basis of the many-body approach are given in Ap- 
pendix|X] We restrict our analysis in this section to the ideal 
limit of Ramsey interferometry, in which the gas is weakly 
interacting during the evolution period. 



1. Atom-molecule interference of a single pair of atoms 

Motivated by the approaches of Refs. I32ll33ll . we first con- 
sider a pair of **^Rb atoms in a spherically symmetric trap, 
exposed to a spatially homogeneous magnetic field, whose 
strength is varied in time, according to the sequence of pulses 
in Fig. [2 The degrees of freedom of the centre of mass and 
relative motions of the atom pair are then exactly decoupled, 
and the dynamics of the relative coordinate r is determined by 
a Hamiltonian of the form: 

H2B = V2 + y,,,p(r) + V{B, r). (2) 

m 

Here VtrapC?") is the potential of the atom trap, m is the atomic 
mass, and V{B, r) is the magnetic field dependent binary in- 
teraction potential. In general, H2B depends on all hyper- 
fine states of the atoms, which are strongly co upled by the 
Feshbach resonance state 0ies('") (cf., e.g., Refs. I20ll32ll and 
Fig. For simplicity, we shall suppress, throughout this pa- 
per, all indices related to the multichannel nature of the bi- 
nary interactions. All our results and derivations are, however, 
quite general and can be applied to zero energy resonances 
with any number of relevant scattering channels. 

The analogue of the condensate mode is realised by the 
lowest energetic state of two atoms above the dissociation 
threshold of the pair interaction potential V{B, r), whose wave 
function 0o(r) thus obeys the stationary Schrodinger equation 
H2B<l>o{f) = Eo(po{r), with the two-body Hamiltonian evalu- 
ated at the initial magnetic field strength of the pulse sequence. 
The amplitude for the probability to detect the atoms in the 
state 00 at the final time ffinai of the pulse sequence thus pro- 
vides the analogue of the amplitude for the remnant conden- 
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sate fraction, and is determined by: 



(3) 



Here t/2B(ffinai, fo) is the time evolution operator, which satis- 
fies the Schrodinger equation 



in — U2B{t,tQ) = H2B{t)U2B(t,to)- 

at 



(4) 



Separating the pulse sequence in Fig. [2 into its different 
components, L'2B(ffinah fo) can be factorised into the evolution 
operators f/#i(fi,fo), C/#2(f final, fa) and f/evoive(fevoive), associ- 
ated with the first and second magnetic field pulse and with 
the evolution period of the sequence, respectively. Here t\ is 
the time immediately after the first magnetic field pulse, f2 is 
the time at the beginning of the second magnetic field pulse, 
and t2-t\ - fevoive is the duration of the evolution period. The 
factorisation of t/2B(ffinah h) yields: 

t/2B(ffinah fo) = t^#2(ffinah ^2) t/evolve(fevolve) (f 1 , fo)- (5) 

During the evolution period the magnetic field strength is con- 
stant and f/evoive( fevoive) depends only on the duration fevoive- It 
can therefore be written in terms of the stationary energy states 
of H2B at the magnetic field strength Bevoive (cf. Fig.Q. This 
leads to the spectral decomposition: 



evolve(fevolve) ^ ^ \4^\' )^ 



"'evolvc/S^^evolve 



(6) 



i'=-i 



Here we have labelled the two-body energy states (^™'™(r) 
by their vibrational quantum numbers, starting from the high- 
est excited vibrational molecular bound state (;A™'™(r) = 
^evoive(^)^ whose energy fibC^evoive) = ^b™'™ depicted in 
Fig. 13 to the quasi continuum of excited states (^™'™(r) 
(v = 0, 1,2, . . .) with energies above the dissociation thresh- 
old of y(fievoive, '")■ As the initial state, as well as the trap 
potential, are spherically symmetric, we have left out all en- 
ergy states corresponding to partial waves with an orbital an- 
gular momentum different from I - 0. We have also neglected 
all vibrational molecular bound states with energies below 
/ibC-Sevoive) because these states are energetically inaccessible. 

Inserting Eq. (|6j into Eq. Q, the transition amplitude in 
Eq. Q can be represented in terms of a diffusion contribution 
£^26 (fevoive), which is slowly varying in fevoive and stems from 
the quasi continuum states, along with a comparatively rapidly 
varying oscillatory amplitude A2B (fevoive) associated with the 
spectral component of the molecular bound state: 



T2B{4>a 0o) - £'2B(fevolve) + A2b( fevoive )• 



(7) 



Both amplitudes can be expressed in terms of the amplitudes 

/ jcvolve 



M =<0r°'"if^#i(fi,fo)i0()), 



^ 2BVVv. 

T*l{<Po <- 0r°'") =<0o|C/#2(ffinal,f2)l0r'"> 



(8) 
(9) 



associated with the probabilities for transitions between (po{r) 
and the vibrational states (^™'™(r) at the magnetic field 



strength Bevoive during the first and second magnetic field 
pulse. This gives the diffusion part to be 



evolve\ —iE^^ 



£'2B(fevolve) - ^ ^^^(^O 

r=0 

T^#l / J evolve , ( \ 

while the oscillatory amplitude is given by 

A 2b( fevoive) 



^evolve 



(10) 



xr2«(C 



evolve-jg-i£™'"fcvoi.c/S 



(11) 



The probability Po,o(fevoive) = 1^26(00 <- h)\ for the two 
atoms to be detected in the state cpoir) at the final time of the 
pulse sequence provides the two-body analogue of the rem- 
nant condensate fraction as a function of fevoive- Neglecting the 
weak dependence of the diffusion term D2B (fevoive) on fevoive, 
^o,o(fevoive) takes the well known form of an interference sig- 
nal: 

f0,0(fevolve)=|£)2Bl^ + |A2Bl^ 

+ 2 |D2bI -y/?^V^sin(Wefevolve + A^). (12) 



Here 



O.b 



b.O 



It#1 / -^evolve 
|7^2B('^b 



jevolve\|- 
Pb >\ 



(13) 
(14) 



are the probabilities for molecular formation and dissociation 
during the first and second magnetic field pulse, respectively, 
and 



(15) 



is the angular frequency of the fringes. The total phase shift 
A(p involves the phases of the complex amplitudes in Eqs. @ 
and (|5Jl, and of D2b- 

The signal in Eq. (I12> has maxima and minima at the evo- 



lution times t™': and f™'". , which fulfil the relations 

evolve evolve' 



sin(c.efeToWe + 



A^) =1, 
Sin(^efe™"lve+A^) = -1, 



(16) 
(17) 



respectively. The visibility VRamsey of the Ramsey fringes pro- 
duced by a single pair of atoms is thus given by: 



V, 



^O.O(fevnlve) 



Ramsey 



^0.0(Co"lve) 



^o,o(f; 



evolve ^ 



evolve-^ 



= 2 /p#i P#2 IV 



.avg 
0,0- 



(18) 



Here Pj^^g - |D2bP + \A.2b^ is the average signal obtained for 
evolution times fevoive, at which the oscillatory part on the right 
hand side of Eq. (I12> vanishes. 

These derivations based on two-body physics reveal that 
the emergence of the Ramsey fringes in Refs. 0,01 does not 
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crucially rely upon the presence of many atoms in a Bose- 
Einstein condensate. The same kind of interferometry could 
be realised, for instance, also with two atoms in a tight atom 
trap of an optical lattice site. These atoms could be bosons 
or non-identical fermions, which both can interact with each 
other via s waves. The visibility of the fringe pattern is de- 
termined by the probability of molecular production in 
the first magnetic field pulse and by the probability P^-y of 
their reconversion into atom pairs in the lowest energetic quasi 
continuum state (po{r), during the second magnetic field pulse. 
The angular frequency of the Ramsey fringes is solely deter- 
mined by the binding energy Ef^™^^^. 



2. Continuum limit 

To extend the two-body considerations to an effective de- 
scription of a homogeneous gas of atoms, we shall first de- 
scribe the limiting process as we move from the discrete trap 
states to a continuum. We thus assume that the two atoms 
are confined to a large volume 'V, whose length scales greatly 
exceed the range of the inter-atomic forces as well as the spa- 

The products 



tial extent of the bound molecular state 0^*°'^'^. 
'VPq'j, and 'VP* q are then virtually independent of 'V, and the 
coupling between the centre of mass and the relative coordi- 
nates of the atom pair is negligible. In the continuum limit, the 
amplitudes r*^(0™'™ f- 0o) and T*l((f>o <- 0^™'^*=) contain 
the two-body time evolution operator and the bound molecular 
state <p^°^'"' in free space, as well as the limit of the lowest en- 
ergetic quasi continuum state \<po} to a zero momentum plane 
wave |0) of the relative motion of the atoms. This implies the 
formal replacement: 



I0o) ^ |0> 



(19) 



This choice of the initial state properly reflects the coherence 
properties of a Bose-Einstein condensate l34ll . 

In contrast to molecular production or dissociation, the 
transfer between quasi continuum states involves a spatial 
range on the order of magnitude of the complete sample vol- 
ume "V. This prevents us from performing the continuum limit 
directly in the elastic amplitude Tr^iipa 0o) of Eq. 0. We 
consider instead the continuum limit of the function 

d 

Kt,to) = 'Vin-{cf>o\U2B(t,to)\cf>o) 
at 

^ {2nn)\0\V{B{t))U2j,{t, fo)|O>0(f - fo), (20) 

which is related to the elastic transition amplitude through the 
following: 



= 1- — 



dth{t,to). 



(21) 



Here 6(t - to) denotes the step function, which yields unity for 
f > f and zero elsewhere. The first term on the right hand side 
of Eq. J21l i can to some extent be interpreted as the part of the 
amplitude associated with the events in which both atoms do 



not scatter, while the second term, proportional to the inverse 
of the volume, accounts for all possible collision events during 
the pulse sequence. 



3. Probabilistic extension of the description of a single pair of 
atoms to the description of a homogeneous gas 

A dilute gas usually occupies a volume 'V much larger than 
the volume needed for a single pair of atoms to collide and 
then separate to a distance outside the range of their binary 
interaction. The volume of the gas can thus be divided into re- 
gions that are small enough to contain at most two atoms but 
still sufficiently large to exceed the range of the inter-atomic 
forces by a large amount. This coarse graining procedure pro- 
vides the physical justification for rescaling two-body transi- 
tion amplitudes in order to describe a homogeneous gas of 
atoms with a density n. The correct scaling factor between the 
effective volume of two colliding atoms and the entire volume 

of the gas can be determined from classical probability the- 
ory. When applied to the molecular conversion efficiency of 
the first magnetic field pulse, for instance, classical probabil- 
ity theory predicts the number 2A^*' of atoms, associated to 
molecules, relative to the total number = of atoms in 
the gas to be: 



N N 



2 N{N - 1) (Infif 



K<^bic/#i(fi,fo)io>r 



,#1 

0,b- 



(22) 



Here 2/N is the number of atomic constituents of a sin- 
gle diatomic molecule relative to the total number of atoms, 
A^(A^ - l)/2 X N- /2is the number of pairs of atoms in the gas, 
and the remaining factor (27:hf K0b|C/#i(fi, fo)|0>P is the 
probability for the association of a single pair of atoms in the 
continuum limit. The probabilistic extension of the descrip- 
tion of a single pair of atoms to the description of a homoge- 
neous gas thus consists in identifying the effective volume of 
a colliding pair of atoms by the inverse of the atomic density. 
This identification is quite intuitive, as 1/n is exactly the vol- 
ume, in which, on average, a single atom can be found in a 
homogeneous gas. We thus replace the two-body amplitude 
on the right hand sides of Eq. i2H by 



2B 



(00 <- 0o) = 1 



■ n— I 



dthitJo). 



(23) 



The probabilistic approach implies, in particular, that we 
should multiply the transition probabilities for molecular pro- 
duction or dissociation in Eqs. (I13> and MA\ by the number 
of atoms = ifV. A similar procedure was previously ap- 
plied by Mies et al. (cf. Refs. I20I l32il '). leading to a re- 
markable agreement between theory and the Feshbach res- 
onance crossing experiments of Ref. lllll . From its deriva- 
tion given in this paper, however, it becomes apparent that 
the probabilistic extension of the two-body physics to the de- 
scription of a gas is not restricted to Bose-Einstein conden- 
sates but could be applied also to the interacting pairs of 
atoms in dilute two component Fermi gases, boson-fermion 
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mixtures or non-condensed Bose gases. A comparison be- 
tween Eq. i23\ and the dynamic equation j28> for the atomic 
mean field of the genuine many-body theory of Refs. 
impHes that the applicabihty of the probabiHstic extension 
of the two-body physics is restricted to time scales ffinai - 
fo « »z/[2;ja(fievoive)nc™''"]> where is the mean con- 

densate density during the evolution period, and 27rS//Yevoive = 
m/[2/!fl(Bevoive)«r°'^''] the oscillation period associated with 
the chemical potential. 




Ql • 1 . 1 • 1 • 1 • 1 . 1 

10 15 20 25 30 35 40 

^evolve ^^^^ 

FIG. 4: Number of remnant condensate and burst atoms in depen- 
dence on ?evolvei determined from the probabilistic extension of the 
two-body dynamics, evaluated at the mean experimental density 0| 
of 3.16 X 10'^ cm"'. All the other physical parameters and symbols 
are the same as those in Fig.|2| 



Figure |3 shows the interference fringes predicted by the 
probabilistic extension of the two-body dynamics for the ex- 
perimental situation described in Fig. 6 of Ref. 0|. The com- 
parison between Fig.|3]of this paper and Fig. 6 of Ref. |0] on 
one hand and Fig.l^on the other hand reveals that the two- 
body approach overestimates molecular production efficien- 
cies and the fringe visibilities in a Bose-Einstein condensate. 
The overall qualitative features, as well as the order of mag- 
nitude of all components, however, correctly reflect the pre- 
dictions of full many-body simulations as well as the experi- 
mental observations. The predictive power of such a simple 
approach based on classical probability theory is even more 
remarkable in view of the fact that the widely used two level 
mean field approach of Refs. lll[ll[ll[lllll[ll[ll fails 
completely to describe the fringe pattern. In fact, this version 
of a many-body mean field theory not only rules out the pro- 
duction of burst atoms I20I1 but it also predicts, when applied 
to the physical situation of Fig. 6 in Ref. 0], that the num- 
ber of condensate atoms should remain virtually unchanged, 
throughout the entire pulse sequence. This failure of the two 
level mean field approach is due to its incomplete description 
of the binary collision physics \20]. 

Equation ( I18> implies that the visibility of the interference 



fringes in the remnant Bose-Einstein condensate is given by 

VRamsey=2«^^P«P*yP^,;^ (24) 

where P"^^ can be obtained from the transition amplitude in 
Eq. (123 > and thus converges to unity in the limit of zero den- 
sity. The product 'V is independent of the density 

n and of the volume Consequently, the visibility VRamsey 
varies linearly with the density in the limit « — > 0. The prob- 
abilistic approach recovers the angular frequency of the Ram- 
sey fringes of the associated two-body problem in Eq. (I15> . 

4. Molecular production efficiency 

The visibility of the Ramsey fringes in Eq. J24> strongly 
depends on the molecular production efficiency Py'j^, during 
the first magnetic field pulse of the sequence in Fig. ^ In 
the course of the molecular association, pairs of atoms sep- 
arated by the average inter-atomic distance of the gas are 
transferred into a highly excited but comparatively tightly 
bound state, with a mean internuclear distance on the order 
of (r) = a(Bevoive)/2. These length scales typically differ by 
several orders of magnitude, and special techniques are re- 
quired to efficiently overcome the problem of the large spatial 
separation of atoms in a dilute gas. 

The technique of magnetic field pulses, applied in Refs. 
01, uses two steps to achieve molecular association. First, the 
magnetic field strength is swept close to the zero energy res- 
onance, so that the scattering length fl(Brnin) becomes compa- 
rable in magnitude to the mean inter-atomic distance of the 
gas. From the intermediate length scale set by the 
atom pairs are then transferred to their mean internuclear dis- 
tance of (r) = fl(Bevoive)/2 in the final upward ramp of the 
first magnetic field pulse. This physical picture suggests that 
an efficient molecular association requires a favourable bal- 
ance between the mean inter-atomic distance of the gas and 
the scattering lengths fl(Bmin) and a(Bevoive)- A similar bal- 
ance of length scales is also required to achieve favourable 
Franck-Condon overlaps between initial and final molecu- 
lar states in the general context of molecular spectroscopy 
|3^ . Our analytic studies of Appendix 151 indeed, indicate 
that the mechanism of molecular production with magnetic 
field pulses exhibits analogies to the experimental technique 
of two colour photo-association (see Ref. IjtIi for a descrip- 
tion of the technique), in which free pairs of atoms are trans- 
ferred into a comparatively tightly bound molecule via Raman 
transitions through a highly excited intermediate state. While 
the photo-association approach can populate a single deeply 
bound target level, by appropriately tuning the laser frequen- 
cies, the molecular production via magnetic field pulses in- 
volves a range of energies that overlap with the continuum 
part of the two-body spectrum. This reduces the efficiency of 
the association of molecules in the highest excited diatomic 
vibrational bound state and leads to the observed production 
of burst atoms 0] . 

These considerations are confirmed by the results on the 
molecular production efficiency in the first pulse of the mag- 
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B . [G] 

mm ^ J 

FIG. 5: The molecular production efficiency P*^^^ as a function of the 
minimum magnetic field strength 6„ii„, during the first pulse of the 
sequence of Fig.Q for a fixed evolution field strength of B„oUc =160 
G. The figure compares molecular production efficiencies of the re- 
alistic pulse of Fig. Q (squares) with predictions associated with an 
idealised pulse sequence (cf. Appendix IaI. which just contains sud- 
den changes of the magnetic field strength. These predictions involve 
calculations using a realistic low energy binary interaction potential 
15, 20] (circles and squares) and the contact potential approx- 
imation of Appendix |B| (solid curve). The duration of the idealised 
first magnetic field pulse of the sequence was chosen to be 36.2 /js. 



netic field sequence presented in Fig.|5] For the case of an ide- 
alised pulse sequence, consisting of sudden jumps of the mag- 
netic field strength, we observe an optimal value of B^in of 
around 156.7 G, where a favourable balance between the three 
inter-particle length scales involved in this process is reached. 
The analogous calculation performed for the realistic pulse 
sequence determines the optimal value of B^m to be smaller 
than 155.7 G, in agreement with the pulse sequence actually 
employed in the experiments of Refs. |7, 8]. Although both 
the realistic and idealised pulse sequence calculations predict 
the same order of magnitude for the molecular production ef- 
ficiency, they significantly differ in their precise dependences 
on the magnetic field Bmin- This reveals the sensitivity of the 
molecular component produced, and, as a consequence, of the 
observed fringe visibility, to the precise shape of the magnetic 
field pulse used. Figure|5]also shows that the analytic estimate 
of the molecular production efficiency in the contact potential 
approximation for the idealised pulse sequence (see Appendix 
|E]for details) recovers both the magnitude and the essential 
features of the predictions based on a realistic low energy bi- 
nary interaction potential. For the values of Bmin smaller than 
156 G we are outside the domain of the analytical approxi- 
mation, as explained in Appendix 151 For the idealised pulse 
sequence we observe effects of quantum interference in the 
inelastic transition amplitude resulting in the oscillations of 
the transition probability for B^in greater than 157 G. These 
interference phenomena in inelastic transition amplitudes are 
sometimes referred to as Stiickelberg oscillations (see, 
e.g., Ref. for similar interferences in three-body recom- 



bination rate constants). 



5. Many-body corrections to the frequency and visibility of the 
Ramsey fringes in the remnant Bose-Einstein condensate 

The main differences between the effective two-body and 
the genuine many-body descriptions of dilute Bose-Einstein 
condensates reside in the treatment of the coherent nature of 
the gas and of the non-linear dynamics of the atomic mean 
field. In Appendix |A] we provide a detailed derivation of the 
many-body corrections to the frequency and visibility of the 
Ramsey fringes in the remnant Bose-Einstein condensate, un- 
der the assumption that the gas is weakly interacting during 
the evolution period of the pulse sequence. In this section we 
give a physical explanation of their origin. For simplicity, we 
shall consider all quantities in the homogeneous gas limit. 

The coherent dynamics of the atomic mean field ^(t), dur- 
ing the evolution period, is largely described by the stationary 
limit ^{t) oc («^™'^'=)'/2exp(-;jUevoivef/S), up to a global phase 
shift, where «^™'™ is the mean condensate density and 

A^evolve = 47rft2fl(B^^„,^^)„-°lve/,„ (25) 

is the chemical potential. The chemical potential can thus be 
largely interpreted as the amount of energy gained by each 
atom due to the presence of the surrounding atoms. We thus 
expect that the finite energy of 2yUevoive in a binary collision 
leads to an upward shift of the two-body angular fringe fre- 
quency in Eq. (I15> . which yields: 

We = (|£b™'"| + 2/«evolve) Ifl. (26) 

This is precisely the estimate that we obtain from the system- 
atic treatment given in AppendixIXl 

We note, however, that the Gross-Pitaevskii approxima- 
tion of the chemical potential in Eq. M5\ becomes singu- 
lar when the magnetic field strength approaches the zero en- 
ergy resonance. This singularity is an artifact of the assump- 
tion of weak interactions, i.e. n^*°'^''a^(Bevoive) ^ 1, during 
the evolution period. The Gross-Pitaevskii approximation of 
Eq. (I25> may be systematically improved using the self con- 
sistent treatment of Ref. fl9ll . which determines the chemical 
potential in terms of the exact stationary solution of the gen- 
uine many-body dynamic equation (I28> for the atomic mean 
field E3l. This stationary approach leads to a chemical po- 
tential, which properly accounts for the energy dependence of 
the binary transition amplitudes lITgil and leads to the physi- 
cally expected smooth behaviour of yUevoive when the magnetic 
field strength is tuned in the close vicinity of the zero energy 
resonance. 

Our studies of SectionH^also indicate that even in the ab- 
sence of deeply inelastic loss processes, like, e.g., spin relax- 
ation, the experimental 0, [3l evolution times fevoive are too 
short for the gas to equilibrate. In particular, fevoive is typi- 
cally much smaller than the time scale ti / (2/ievoive) set by the 
many-body frequency shift in Eq. ( I26> . As a consequence, the 
experiments |7j,Q] could not fully resolve the frequency shift. 
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The derivations in Appendix IaI reveal that the many -body 
corrections to the visibility in Eq. (I24> are mainly due to the 
non-linear dynamics of the atomic mean field. The final result 



1 



Ramsey 



■exp(-4n-F[pg'jt:o]"') 
.exp(-4„n/[p«P*5,]'^') 



(27) 



largely recovers the dependences of Eq. J24l i on the physical 
parameters n, ^Pq\ and 'VP^ so that both formulae asymp- 
totically coincide in the limit n 0. The non-linear nature 
of the mean field dynamics, however, leads to a reduction of 
the fringe visibility in comparison with the probabilistic ex- 
tension of the linear two-body time evolution, when the den- 
sity is increased. A similar behaviour can be observed in the 
comparison between genuine many-body predictions and the 
extended two-body Landau-Zener approach of Ref. ll32il for 
the description of the molecular production in Feshbach reso- 
nance crossin g ex periments with linear ramps of the magnetic 
field strength Boll. 



III. DECOHERENCE OF THE RAMSEY FRINGES 

In this section we consider in detail the regime of magnetic 
field strengths in the close vicinity of the zero energy reso- 
nance, in which decoherence phenomena diminish the Ram- 
sey fringes in the remnant Bose-Einstein condensate. This 
regime of low molecular vibrational frequencies is strongly 
influenced by the non-linear dynamics of the atomic mean 
field. Our numerical simulations are therefore based on the 
many-body approach of Refs. 0, fioll that includes both a 
non-perturbative treatment of the binary physics and the non- 
linearities to all orders in the atomic mean field. 



A. Comparison between experimental Ramsey fringes in the 
remnant Bose-Einstein condensate and many-body simulations 

In the first order microscopic quantum dynamics approach 
of Refs. llgl fioll the evolution of all components of the gas is 
determined by a single non-linear Schrodinger equation for 
the atomic mean field, i.e. 



at 



v^ + y„,p(x) 

2m 



*l'(x, f) 



dT^'^(x,T) — h{t,T). (28) 
OT 



The inter-atomic collisions are included fully in this approach 
in the time dependent coupling function of Eq. i2Q\ . The dy- 
namics of the atomic mean field determines not only the con- 
densate density «c(x, t) = |' y(x , f)P but also the molecular and 
burst fractions of the gas llOll . The simulations of the final 
number of atoms of all components of the gas in Fig.|2lcorre- 
spond to the experimental situation of Fig. 6 in Ref. |7], and 
were obtained from a sequence of full solutions of Eq. j28> 
with variable evolution times. 



In the course of our studies, we have solved Eq. ( I28> for a 
variety of experimental |@1 evolution times fevoive and mag- 
netic field strengths Bevoive- Our simulations fully account 
for realistic sequences of magnetic field pulses 1251 . shown 
schematically in Fig. [3 and for the inhomogeneous nature of 
the Bose-Einstein condensate in a cylindrical atom trap with 
the experimental radial and axial frequencies of Vradiai = 17.5 
Hz and Vaxiai - 6.8 Hz, respectively. The initial state of 
our simulations is the ground state of the stationary Gross- 
Pitaevskii equation, corresponding to a '^^Rb Bose-Einstein 
condensate with 16000 atoms j^l at the initial magnetic field 
strength of Z? = 162.2 G of the pulse sequence in Fig.^ 
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FIG. 6: Comparison of the experimental IM l2^ number of atoms in 
the remnant Bose -Einstein condensate as a function of fcvoive (cir- 
cles) with the predictions of Eq. <28> (squares) at the end of each 
magnetic field pulse sequence. The dashed and sohd curves indicate 
the fits of Eq. <29> 01 to the experimental and theoretical data, re- 
spectively. In the course of the experiments the magnetic field 
strength ficvoivc during the evolution period of the pulse sequence was 
varied between 156.076 and 161.771 G. The fringe patterns cover 
a wide range of evolution times, from about lOyus at 161.771 G to 
about 1 ms at 156.076 G. 



Figure |6] shows comparisons between the experimentally 
observed dependence on fevoive of the remnant condensate and 
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theoretical predictions of Eq. ( I28> . at the end of each magnetic 
field pulse sequence. Both the experimental data and the the- 
oretical predictions exhibit damped oscillations, whose decay 
becomes particularly marked for sequences with longer evo- 
lution times on the order of 1 ms, at magnetic field strengths 
of Bevoive S 157 G, i.e. in the close vicinity of the zero energy 
resonance at Bo - 155.041 G jlsl). The theoretical predictions 
have the same systematic trends as the experimental fringe 
patterns with frequencies between 4 kHz at 156.076 G and 
about 1 MHz at 161.771 G. All of the theoretical curves over- 
estimate the measured numbers of remnant condensate atoms. 
These deviations may, at least in principle, be due to fast in- 
elastic loss phenomena such as spin relaxation in **^Rb col- 
lisions, whose rates are enhanced in the vicinity of the zero 
energy resonance Q- We shall show in Subsection 1111 CI 
however, that spin relaxation is unlikely to play a dominant 
role in these deviations. Our results of Subsection lll CI rather 
indicate that the remnant condensate fractions in the interfer- 
ometric experiments are very sensitive to the precise shape of 
the magnetic field pulse sequence. We therefore believe that a 
mismatch between the real experimental pulse sequence and 
the sequence used in our simulations is a more likely reason 
for the deviations, in particular, for magnetic field strengths 
^evolve ^ 158 G, away from the position of the zero energy 
resonance. 



B. Oscillation frequencies of the Ramsey fringes in the 
remnant Bose-Einstein condensate 

In order to extract the fringe frequencies from our simu- 
lations, we have closely followed the procedure suggested in 
Ref. I §S and fitted the theoretically predicted dependence of 
the number of remnant condensate atoms on fevoive to the for- 
mula: 

Arevolve _»r _ n,f 

~^*avg ^'evolve 

+ Ae^^'""'" sin (Wefevolve + A^). (29) 

The angular frequency Wg of the fringes has been interpreted 
in Ref. |8J in terms of a natural oscillation frequency vq, 
which is related to the angular frequency oje through We = 

In - \J3l{2n)Y. The amplitude A of the interference 

fringes as well as the average number A^ave of atoms deter- 
mine the fringe visibility. The remaining fit parameters are 
a loss rate of atoms a, a damping rate /? and a phase shift 
A(/? of the Ramsey fringes. We note that the fit parameters of 
Eq. \29\ can be largely identified from the interference sig- 
nal of a single pair of atoms in Eq. M2\ . provided that the 
weak dependence of the diffusion contribution Z)2b (fevoive) on 
the evolution time fevoive is taken into account. The fit at 
^evolve - 156.076 G in Fig. |5] indicates, however, that the 
functional form of Eq. ( I29> is not entirely compatible with 
the non-linear dynamics of the atomic mean field in Eq. (I28> . 

Figure shows the natural frequency vq in dependence 
on Bevoive obtained from the fits to the experimental and 
theoretical data in Fig. |6] respectively. As pointed out in 
Refs. i7 nl0i.i26i1 . for magnetic field strengths Bevoive > 158 G, 




FIG. 7: Natural frequency vo of the Ramsey fringes as a function 
of the magnetic field strength Bcvoive as obtained through fits of the 
theoretical predictions (squares) and of the experimental I^ E3l num- 
ber of remnant condensate atoms (circles) to Eq. <29> . The solid and 
dashed curves indicate theoretical predictions of the fringe frequency 
that account for the first order many-body contribution 2/ievoivc/'i 
(cf. Subsection IllO of the chemical potential in the self consis- 
tent stationary approach of Ref. fl^ (solid curve) and in the Gross- 
Pitaevskii approximation of Eq. <25> (dashed curve), respectively. 
The dotted dashed curve shows the binding frequency of the highest 
excited vibrational molecular bound state of '*'Rb2, in dependence on 
the magnetic field strength. The inset shows an enlargement of the 
near resonance region of magnetic field strengths Sevoive, in which 
the deviations between fringe frequency and molecular vibrational 
frequency become apparent. Small error bars in the theoretical data 
result from the uncertainties associated with the fitting procedure. 



sufficiently far away from the zero energy resonance at Bq = 
155.041 G, the fringe frequencies closely match the vibra- 
tional frequencies associated with the weakly bound molec- 
ular state 0^™'^''. A comparison with Fig. |3] shows that the 
binding energies in this range of magnetic field strengths are 
sensitive to the long range van der Waals dispersion coefficient 
C(, fioll and even to the magnetic moment dEres/dB of the 
Feshbach resonance level (cf.,e.g., Refs. 12(11321 '). which de- 
termines the linear slope of E\,{B) at magnetic field strengths 
^evolve > 162 G. 

Closer to the zero energy resonance both the experimental 
data 01 as well as the theoretical predictions show a visible 
upward shift of the fringe frequencies from the molecular vi- 
brational frequencies. From the analytic results of our theoret- 
ical studies in Section HI] and from our numerical simulations 
we expect this shift to persist for all values of Bevoive but to de- 
crease in its relative magnitude with increasing magnetic field 
strengths. The four experimental points closest to the zero 
energy resonance, to which a shift has been attributed in the 
experimental reference 1 8], agree in their magnitudes with our 
theoretical predictions (see inset of Fig.0. For magnetic field 
strengths Bevoive ^ 156.84 G the shift was neglected in Ref. 0] 
and the associated fringe frequencies were used to adjust the 
theoretical curve of the binding frequencies in Figs.|31and0 
This, in turn, also determines the parameters Bq, AB and Obe 
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of Eq. Q. As we have used these parameters as an input to 
the coupling function h(t, r) of Eq. j28> , and the frequency 
shift persists for all magnetic field strengths Bevoive in our ap- 
proach, we expect our simulations to deviate slightly from the 
experimental findings simply due to the slight uncertainty in 
the Feshbach resonance parameters. 



C. Two-body decay of the Bose-Einstein condensate 



Despite the fact that the fit parameter a of Eq. ( I29t is de- 
termined with a much larger uncertainty than the frequency 
vo, we may interpret it in terms of two-body decay processes 
during the evolution period. Two-body loss phenomena in a 
Bose-Einstein condensate are quite generally described by the 
rate equation 



A^c(0 = - Y<«c(f))A^c(0. 



(30) 



Here K2 is the two-body loss rate constant, {n^it)) is the aver- 
age condensate density, and A^c(0 is the number of condensate 
atoms. We note that Eq. j30> is not restricted to the descrip- 
tion of spin relaxation events but refers to all two-body loss 
phenomena which occur on time scales much shorter than 
those set by the bulk motion of the gas. We shall restrict 
our considerations on the two-body decay to magnetic field 
strengths Bevoive < 157 G in the close vicinity of the zero en- 
ergy resonance, where we find that the pulse sequences are 
sufficiently long to provide reliable values for a with a sys- 
tematic trend as a function of Bevoive- In this region the Ram- 
sey fringes of Fig. |6l also exhibit a strong decoherence, and 
diatomic molecules in the highest excited vibrational bound 
state can not be identified as a separate entity because their 
wave functions would be more extended than the mean inter- 
atomic distance of the gas ITlll . In order to establish the rela- 
tionship between a and K2, we solve Eq. ( I30> formally by di- 
viding both sides by N^it) and then integrating the logarithmic 
derivative from the time fi, immediately after the first mag- 
netic field pulse of Fig.Q] to the time f2, immediately after the 
evolution period. This yields: 



A^c(f2) = A^c(fi)exp - 



r dt{n, 
Jt, 



(t)} 



(31) 



If we assume that during the evolution period fevoive - h - 
ti the loss of condensate atoms is small in comparison with 
Nc(t\), we can expand the right hand side of Eq. (13 U to first 
order in fevoive- This leads to the following expression: 

A^c(f2) ^ A^c(fl) - YA^c(fl)<«c(fl))fevolve. (32) 

A comparison between Eqs. \29l and J32> then allows us to 
identify the parameters A^c(fi) - A^avg, A^c(r2) = A?™'™ and 



K2 = 



2a 



N^wgiriciti))' 



(33) 



provided that the second magnetic field pulse does not sig- 
nificantly change the number of condensate atoms. Our nu- 
merical simulations show that this assumption is justified for 



longer pulse sequences with evolution fields devolve < 157 G 
(cf. Fig.|5), while for the faster sequences it may be violated. 
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FIG. 8: The two-body loss rate constant K2 determined from the fit 
parameter a of Eq. <29> through Eq. <33> as a function of the mag- 
netic field strength Sevoivc during the evolution period of the pulse 
sequence. The circles indicate the experimental values (igj, while the 
squares correspond to the theoretical predictions of Eq. <28> depicted 
in Fig.|6| The error bars indicate the statistical errors of the fits. The 
solid curve indicates coupled channels calculations iiHl of the spin 
relaxation rate constant for a pair of '^'Rb atoms in the limit of zero 
collision energy. The rate constant K2 is given on a logarithmic scale. 



Figure|S]provides an overview of the magnitudes of the two- 
body loss rate constants K2 obtained through Eq. j33t from 
the fits to the experimental data |3| (circles) and to the theo- 
retical predictions (squares) of Fig.|6l The rate constants K2 
are shown as a function of the magnetic field strength Bevoive 
during the evolution period. For comparison, the solid curve 
indicates results of a coupled channels calculation L41.1 of the 
spin relaxation rate constant. The comparatively large size 
of about 10"^ cm-'/s of both the experimental (sl] and theo- 
retical rate constants strongly indicate that the loss of con- 
densate atoms is mainly due to the energy transfer from the 
first magnetic field pulse. This drives the initially weakly in- 
teracting Bose-Einstein condensate into a strongly correlated 
non-equilibrium state. Deeply inelastic spin relaxation loss 
phenomena are not included in our simulations of Eq. ( I28t . 
The theory thus predicts that the atoms are transferred from 
the condensate into correlated pairs of atoms in excited states. 
The coupled channels calculations Il41il . indeed, confirm that 
the rate constants for spin relaxation (sohd curve in Fig.|8} are 
less than 10"" cmVs in the relevant range of magnetic field 
strengths, i.e. from 156 to 158.5 G, and thus do not explain 
the size of the observed losses. Our interpretation is corrobo- 
rated by the observation of Ref. 1 8] that over the whole range 
of magnetic field strengths Bevoive the loss rate a is weakly 
dependent on Bevoive and consistent with the number loss in a 
single magnetic field pulse 142] . 
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D. 



Visibility of tlie Ramsey fringes in the remnant 
Bose-Einstein condensate 
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FIG. 9: Visibility of the Ramsey fringes in the remnant Bose-Einstein 
condensate determined from the experimental data (circles) and from 
the solutions of Eq. <28> (squares), as a function of the magnetic field 
strength Bcvoivc during the evolution period of the pulse sequence. 
We have also plotted the estimate of the fringe visibility given by 
Eq. )27> . evaluated for the magnetic field pulses of the realistic se- 
quence (diamonds) and of an idealised pulse sequence with sudden 
changes of the magnetic field strength (cf. Appendix IaI (triangles). 
For the idealised pulse sequence the durations of the first pulse and 
the second pulse were chosen to be 36.2 /is and 23.9 yus, respectively. 
Equation <27t has then been evaluated for n = 3.07 x 10'^ cm"^, 
which corresponds to the initial mean density of the condensate in 
the experiments of Ref. |.80. The dashed curve shows the predictions 
of the contact potential approximation to the probabilities P^^^^ and 



Figure |5] shows the fringe visibility in the remnant Bose- 
Einstein condensate. This is determined as a function of the 
magnetic field strength Bevoive in the evolution period from 
the experimental data and from the full solutions of Eq. j28t . 
The visibilities were obtained from the fits of Eq. j29t to the 
experimental and numerical data of Fig.|6lusing the relation 



VRamsey - A/N; 



avg- 



(34) 



The theoretical results follow the rough size and mono tonic 
trend of the experimental visibilities. The increase of VRamsey 
with decreasing Bevoive is due to the enhanced depletion of the 
average number A^avg of remnant condensate atoms at mag- 
netic field strengths close to the zero energy resonance. This 
trend can also be observed in the visibilities estimated from 
Eq. i27\ . These estimates were obtained from the exact tran- 
sition probabilities P^^^ and P* ^ for both the pulses of the real- 
istic sequence (as depicted in Figure[3 and their counterparts 
in an idealised pulse sequence of sudden changes of the mag- 
netic field strength (cf. AppendixlAl. FigureOalso shows an 
estimate of VRamsey based on Eq. i21l and the contact poten- 
tial approximation (cf. Appendix |b}. The contact potential 



approximation accounts for the inter-atomic interactions just 
in terms of a single parameter, the scattering length a{B). The 
large difference between the contact potential approximation 
and the experimental data as well as the exact theoretical data 
reveals the sensitivity of the interferometric experiments of 
Refs. 0, @] to the nature of inter-atomic interactions beyond 
universal considerations, which are based solely on the scat- 
tering length. This sensitivity of the Ramsey interferometry 
with atoms and molecules Q |8|] to non-universal properties 
of the binary collision physics provides a crucial probe of the 
microscopic physics through the bulk properties of the gas. 
The information it provides is largely inaccessible to exper- 
imental studies involving linear ramps of the magnetic field 
strength across zero energy resonances fl^l20ll . 



IV. CONCLUSIONS 

Our theoretical studies in this paper reveal the remark- 
able range of physical regimes that are probed by the atom- 
molecule coherence experiments of Refs. 0, These 
regimes include the Ramsey interferometer limit, at evolution 
fields sufficiently far away from the zero energy resonance for 
the binary and bulk motions to be virtually decoupled. In the 
other limit, i.e. the close vicinity of the zero energy resonance, 
two- and many-body phenomena can not be disentangled dur- 
ing the evolution period. 

Using a probabilistic extension of the dynamics of a single 
pair of atoms to the description of a gas, we have clearly iden- 
tified the different elements of a Ramsey interferometer from 
the experimental OilH sequences of magnetic field pulses. We 
have shown that much of the observations of Refs. 001, in 
particular, the frequency and visibility of the Ramsey fringes, 
can be understood qualitatively and also largely quantitatively 
on the basis of two-body physics. Our studies reveal that the 
probabilistic extension of the pair dynamics is applicable pro- 
vided that the experiments are operated in the interferometer 
mode, i.e. at evolution fields far away from the zero energy 
resonance. These studies were motivated by the approaches 
of Refs. I3^l33ll . Our probabilistic approach, however, is not 
restricted to Bose-Einstein condensates and does not include 
any adjustable parameters. 

We have shown that in the interferometer mode, in con- 
tradiction to the interpretations of some previous theoretical 
studies the gas is virtually unaffected by non- 

equilibrium phenomena during the evolution period of the 
pulse sequence, and the mutually orthogonal molecular and 
free atomic components evolve independently and coherently. 
The Ramsey fringes thus reflect the phase difference of both 
components rather than dynamical exchange of particles be- 
tween them. In fact, such an exchange between the molec- 
ular and free atomic components during the evolution period 
would imply the breakup and reformation of molecules. In the 
absence of any external driving force these would be energet- 
ically impossible. 

In the Ramsey interferometer mode the measurements pro- 
vide detailed insight into the binary collision dynamics be- 
yond universal considerations. These interferometric expert- 
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ments are thus sensitive to the long range van der Waals dis- 
persion coefficient Ce of the binary interactions and even to 
the magnetic moment of the Feshbach resonance level, rather 
than just to the scattering length. Our studies clearly show that 
this sensitivity is observed not only in the fringe frequencies 
but also in their visibilities. The existence of the fringe pattern 
and its dependence on the physical parameters of the system 
is remarkable in its own right with the experiments clearly 
showing how vibration of molecules on a nanometre length 
scale can significantly influence the dynamics of the entire gas 
on length scales as large as several micrometres. This access 
to the microscopic physics can not easily be gained from the 
molecular association techniques based on linear ramps of the 
magnetic field strength across zero energy resonances, which 
are mainly sensitive just to universal scattering observables 
MM- 

Operated in the close vicinity of the zero energy resonance, 
the gas is driven into a strongly correlated non-equilibrium 
state during the evolution period. In this regime the binary 
collision physics and many-body phenomena can not be dis- 
entangled because the wave function of the highest excited 
diatomic vibrational bound state strongly overlaps with the 
mean distance between the atoms in the dilute gas ifTlll . A 
theoretical description of the observations of Ref. (31 thus 
requires genuine many-body considerations, which fully ac- 
count not only for the two-body physics but also for the non- 
linear dynamics of the atomic Bose-Einstein condensate in a 
non-perturbative manner Our dynamical simulations of all 
experiments of Ref. were based on the microscopic quan- 
tum dynamics approach to dilute Bose-Einstein condensates 
of Refs. [Toll and explained all the systematic trends that 
we were able to identify from the measured data. These trends 
involve an upward frequency shift of the interference fringes, 
accompanied by decoherence and an overall decay of the con- 
densate. We have shown that the degree to which these phe- 
nomena influence the fringe pattern can be interpreted as a 
measure of the extent to which molecules in the highest ex- 
cited diatomic vibrational bound state can be identified as a 
separate entity of the gas. Our simulations indicate that even 
in the physical regime close to the zero energy resonance, in 
which the pulse sequences are sufficiently long to reveal a de- 
cay of the condensate, the non-equilibrium dynamics is deter- 
mined mainly by the same physical parameters of the binary 
interaction potential as in the interferometer mode. Deeply 
inelastic loss processes, such as spin relaxation, are likely to 
be negligible. Such processes will, however, affect the parts 
of the experiment after the magnetic field pulse sequence, in- 
volving the ballistic expansion of the gas before imaging. On 
such longer time scales spin relaxation leads to the sponta- 
neous dissociation of the molecules into fast fragments, which 
explains the presence of the component of undetected atoms 
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APPENDIX A: MANY-BODY DYNAMICS OF THE 
BOSE-EINSTEIN CONDENSATE 

In this appendix we derive the formulae for the frequency 
and visibility of the Ramsey fringes in the remnant Bose- 
Einstein condensate [cf. Eqs. (I26t and ( I27H from the first or- 
der microscopic quantum dynamics approach of Refs. Igl fioll . 
We restrict our considerations to the region of magnetic field 
strengths Bevoive, in which the gas is weakly interacting during 
the evolution period of the pulse sequence. In this ideal limit 
of the Ramsey interferometry with atoms and molecules in a 
Bose-Einstein condensate the phenomena associated with the 
non-linear dynamics of the atomic mean field are weak. For 
simplicity, we perform all derivations in the homogeneous gas 
limit. 



1. Hydrodynamic representation 

In the homogeneous gas limit the atomic mean field *P(f) 
depends just on the time variable t. To study the interference 
fringes, it is convenient to represent in terms of its uni- 
form density and phase; 

^PCO = ^fh(fje-'^^^'\ (Al) 

In this representation we then introduce the exponential form 

«c(f) = (A2) 

of the density. The evolution of the exponent y(f) as well as 
the phase (p{t) are determined simultaneously by the real and 
imaginary parts of the dynamic equation for the atomic mean 
field, i.e. by 

which can be derived directly from Eq. j28> in the homoge- 
neous gas limit. 

2. Idealised pulse sequence 

For simplicity, we idealise the pulse sequence in Fig.^by 
replacing the linear ramps of the first and second pulse by 
sudden changes of the magnetic field strength B{t) at times 
fo < ?! < f2 < ^3 = ffinah where fi - t() and - ti are 
the durations of the first and second pulse, respectively, and 
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fevoive = h - h is the evolution time. The two-body time evo- 
lution operator in Eq. (|5} then factorises into the evolution 
operators at constant magnetic field strength: 



£/2B(f3,fo) = U3{t3 - t2)U2(t2 ' hWiiti - to). 



(A4) 



The associated magnetic field strengths Bi = Bmin = ^3 and 
Bi - devolve are shown in Fig. ^ At each magnetic field 
strength B, (/ = 1,2,3) of the pulse sequence, the two-body 
Hamiltonian is stationary, and Ui(t - can be expanded in 
terms of the stationary energy states: 



f/,(f-f,-i) ^MBi)}e 



-iE^,{Bi){t-l,-i)lh 



I 



dp \4>'^;\Bi)}e 



{MBd\ 
-'■^"-'-)/^'<<'(B,)|. (A5) 



Here (pbiB,) is the the highest excited vibrational bound state 
whose binding energy we have denoted by Eb(Bi), and (^p^'(Z?,) 
is the continuum energy state associated with the momentum 



J 



p of the relative motion of a pair of atoms f4-4\. Through- 
out this appendix we choose the continuum wave functions to 
behave at asymptotically large inter-atomic distances r like 



(r) ~ (Inn) 



-3/2 



(A6) 



where cos i? = p ■ r/{pr) is the scattering angle, and f{p, is 
referred to as the scattering amplitude. The energy associated 
with cpp (Bi) is given by p /m. 

To analyse the evolution of the phase if(t) as well as the 
exponent y(f), we shall successively consider each period of 
constant magnetic field strength B, (/ =1,2, 3). The dynamic 
equation jA3> can be integrated from the time f,_i of the previ- 
ous sudden change of the magnetic field strength to the actual 
time f, which thus fulfils f,_i < f < f, . A separation of Eq. iA3\ 
into its real and imaginary parts then yields: 



<f(t) = ^(f,_i) + Irc j J' dt' -j dT ±h.(t' - r) + C(f', f,_i) j , 

1 / r /^co o 

r(f) = r(fi-i) - -Im dt' e^'VC) _ g-2[y(r)+,V(T)] _ ^) _^ 



Ji-i) 



r 



(A7) 
(A8) 



Here -t) = i2jih)^{0\V(Bi)Ui(t' -T)\0}6(t' -t) is the cou- 
pling function of the dynamic equation JA3> restricted to the 
actual period f,_i < t' < ti of constant magnetic field strength 
Bi, and 



tion 



C(f',f,_,) 



r 



dT^\T)^h{f,T) (A9) 
OT 



accounts for two-body correlations produced before the previ- 
ous sudden change of the magnetic field strength at time f,_i . 



3. Iterative solution of tlie non-linear Schrodinger equation 

We shall show that an approximation based on iteration of 
Eqs. ( IA7> and (IA8> is sufficient to recover the essential quali- 
tative phenomena associated with the dynamics of the atomic 
mean field We thus consider the time integrals involv- 
ing the coupling function /i,(f' - t) and the correlation term 
C(t' , ti-i) as a small perturbation in comparison with the ini- 
tial phase ip{ti-i) and the exponent y(f,-i). 

The first order perturbation approximation consists in re- 
placing, on the right hand sides of Eqs. iA7\ and jA8> , 
the phases ip{t') and ^(t) as well as the density «c(t) - 
exp(-2y(T)) by their initial values (f{ti-\) and ndti-i) - 
exp(-2y(f,_i)), respectively. This procedure then shows that 
in the periods of constant magnetic field strength the approxi- 
mate atomic mean field ^'{t) is subject to a linear time evolu- 



"Vit) « expj-^ J dt' //eff(f',f,-i)j*l'(r,-i). (AlO) 

with a complex one dimensional eff'ective Hamiltonian 

H,ff(t', ti-i) = «c(f,-i)/j,(f' - + e'''^^''-''C(f', f,-i) (All) 

This gives the following expression for the density of the 
atomic Bose-Einstein condensate: 



«c(0 = exp 



lm[//eif(f',f,-i)] hc(f,-i). (A12) 



The dynamics of the condensate mean field is driven by the 
two-body correlations in Eq. iA9\ . which leads to the explicit 
time dependence of //eff(f', f,_i). Equations jA9t and ( lAlOt 
can be used to produce an iterative approximate solution of 
the non-linear dynamic equation iA3i for the entire pulse se- 
quence, just on the basis of the two-body coupling function. 
Figure ^| shows a comparison between the dynamics of the 
condensate predicted by the iterative approximation and the 
exact solution of Eq. ( IA3t . 

The first contribution to the sum on the right hand side of 
Eq. (I^nj provides the initially uncorrected part of the ef- 
fective Hamiltonian. Its contribution to the exponent of the 
atomic mean field ^'{t) in Eq. ( lAlOt can be estimated analyti- 
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FIG. 10: Comparison between the dynamics of the atomic conden- 
sate as obtained from the exact numerical solution (dashed curve) of 
the non-linear Schrodinger equation <A3> and its first order iteration 
approximation (solid curve) for the idealised pulse sequence consist- 
ing of three periods with constant magnetic field strength. The dy- 
namics has been determined for a magnetic field strength in the evo- 
lution period of Sevoive = 160 G, an evolution time of fevoive = SOyUS, 
and an initial condensate density of n = 3.16 x 10'^ cm"'. The inset 
shows a comparison between the Ramsey fringes in the remnant con- 
densate as determined from the two different approaches to Eq. {AH . 



cally (cf. Appendix|Bj from Eq. iA5\ to be 

-n^{ti-\)- dT hi(T) =A,(f - f,-i) + Di(t - f,_i) 

^ Jo 

-iMt-ti-i)/n (A13) 

for each period of constant magnetic field strength B,. Here, 
the leading contribution 



fii = n^{ti_x)a{Bi) 

m 



(A14) 



is the mean field chemical potential of the Gross-Pitaevskii 
equation, 

Ai(t - f,_i) = 87r«e(f,-i)fl'(B,) [e-'-^^(«'K'-'->)/f' - 1 j (A15) 

is an oscillating contribution associated with the bound state 
spectral component of Eq. iA5i . and 

Di(t - r,-i) = - 47r«,(f,_i)fl3(B,){2 ^Ji\E^(Bi)\(t-ti^l)/(7Tr^) 
-[l -w{i^Ji\E^,(Bi)\{t - ti^i)/n)]\ (A16) 



describes a diffusion in density and phase of the atomic mean 
field. The latter stems directly from the continuum part of 
Eq. iA5i . Here, we have chosen the complex square root 
V/ = exp(;.7r/4) and the function denotes the complex 
continuation of the error function H3\ : 



w(z) - e -■"erfc(-/z) 



J z- 



X 



(A17) 



The binding energies have been estimated in terms of the 
scattering length fl(B,) by the universal formula E\,{Bi) = 
-n^l[ma\Bi)l 

The integral involving the initial correlations C(f', fo) van- 
ishes [cf. Eq. ( IA9H and Eq. ( IA13t is, therefore, sufficient to 
determine the evolution of *l'(f) during the period to < t < ti 
of the first magnetic field pulse. In our approximate treatment 
the atomic mean field at the the beginning of the evolution 
period of the idealised pulse sequence is thus given by; 



^(fi) 



(A18) 



Moreover, since we have assumed the gas to be weakly in- 
teracting during the evolution period, it can be verified from 
Eqs. JAUI i. JA14I I. ( IA15> and ( IA16t that within this period of 
time t\ < t < t2 the atomic mean field in Eq. ( lAlOt is well 
approximated by the solution 



>P(f) = e 



-ifi2(t-t,)/n 



(A19) 



of the Gross-Pitaevskii equation, with the initial condensate 
wave function *P(fi) given by Eq. JA18> . 

We shall focus in the following on the oscillatory depen- 
dence of the remnant condensate density «c(f3) in Eq. (IA12t 
on the evolution time fevoive = h-ti, which determines the vis- 
ibility and the frequency of the Ramsey fringes. To this end, 
we introduce the spectral decomposition of the two-body time 
evolution operator 1/2(12 - fi) [cf. Eqs. JA4> and ( IA5H in the 
evolution period into the, as yet undetermined, initially corre- 
lated part of the effective Hamiltonian H^sit' , t2) in Eq. ( lAl H . 
evaluated at times t2 < t' < f3 during the final magnetic field 
pulse of the sequence. A simple calculation then shows that 
the associated contribution to the exponent of the final atomic 
mean field ^{tj,) in Eq. ( lAlOt can be separated into high and 
low frequency components: 



dt e2'V('2)C(f, t2) = A(fevolve) + i'Cfevolve). (A20) 



The high frequency part, i.e. 

X V2 >Pb(ri + 0)e'l^'''*^>l'™'-/^ (A21) 

is associated with the spectral component of C/2(f2 - h) in 
the molecular bound state 0b(B2), and involves the molecular 
mean field fl(ill 



+ 0) 



dT^'-(T) — 
J,„ dT 

^^^^^'^tiXMBim 
V2 



<<^b(52)|t/i(fi -t)|0> 
(A22) 



at the beginning of the evolution period. Furthermore, the ma- 
trix element of U^itT, - t2) on the right hand side of Eq. JA21l i 
can be identified with the two-body probability amplitude 
for transitions from \(l)b{B2)) to the zero energy mode |0o) = 
|0) -\/(27rtiy/'V during the second, final magnetic field pulse 
of the sequence. The probability for the dissociation of a 
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molecule in the bound state 0^^°'^'^ into a pair of condensate 
atoms during the second magnetic field pulse is given by P* q 
[cf. Eq. (O]. 

The remaining low frequency component Z)(fevoive) on the 
right hand side of Eq. ( IA20> involves an integral over relative 
momenta p of the amplitude associated with the burst compo- 
nent of the gas fl^, which has evolved until the beginning of 
the evolution period. This diffusion part Z)(fevoive) is thus de- 
termined mainly by the continuum part of 1/2(12 - h)- Its de- 
pendence on the evolution time fevoive is comparatively weak. 
This allows us to identify the frequency as well as the visibil- 
ity of the Ramsey fringes in the remnant condensate density 
from A(fevoive) in Eq- ( IA2U . while Z)(fevoive) contributes to the 
average of the remnant condensate density with respect to the 
evolution time fevoive . as well as to the small decoherence of 
the interference fringes. We shall focus in the following on 
the determination of the contributions to A(/evoive) on the right 
hand side of Eq. JA21> . 



4. Molecular mean field in the evolution period of the pulse 
sequence 

To study the density and phase shift associated with the 
molecular mean field ^'h{t\ + 0) in Eq. ( IA2U . we perform a 
partial integration in Eq. JA22l i. which yields: 



^b(fi + 0) = 



3/2 w 



V2 

Fv 



r dTmB2)\U^(h-Tm^^\T) 
J la 



dT 



-r2«(C°'^^^0o)^'(fo). (A23) 

The second contribution to the sum on the right hand side 
of Eq. ( IA23t involves the probability amplitude for the as- 
sociation of two asymptotically free zero energy ground state 
atoms into a molecular bound state 0b(fi2) = 0^^°'^'^, during 
the first magnetic field pulse of the sequence [cf. Eq. ©]. The 
first contribution, involving the time derivative of the atomic 
mean field, turns out to be negligible. The molecular density 
^evolve _ |v{/|^(fj + 0)1^ then recovers the result we found using 
the binary classical probability theory: 



•"V = NZ 



— P" 
2 



(A24) 



Here A? = «c(fo)'V and A^J 



evolve 
b 



of atoms and the number of molecules in the evolution period, 
respectively, A^^/2 ^ A^(A^ - l)/2 is the number of pairs of 
atoms in the gas, and is the probability for the association 
of a single pair of condensate atoms with the first magnetic 
field pulse of the sequence to a diatomic molecule in the state 



^evolve [^f £q JJ21l]. 



5. Frequency of the Ramsey fringes in the remnant condensate 

In the first order iterative approximation to the non-linear 
Schrodinger equation iA3\ the remnant condensate density 



«j,(f3) = /ic'^'Vfevoive) is determined by Eqs. ( IA12> . ( IA18> . 
^T^ . and iKM to be 



n^'^fevoive) =ne 



2Re[^3 (rj -h. )+D, (t, -ti)] „2Re[A(/evolve)+£'(levolve )] 



,2Re[A|(/,-f„)+D,(r,-/o)] 



(A25) 



where n - nc(fo) is the total density of the gas. Since the dif- 
fusion term Z)(fevoive) is weakly dependent on fevoive, the oscil- 
latory amplitude A(fevoive) determines the frequency and vis- 
ibility of the Ramsey interference fringes. We shall neglect 
the dependence of Z)(fevoive) on fevoive in the following consid- 
erations and represent the two-body transition amplitudes of 
Eqs. (jSji and (|9jl in terms of their associated probabilities and 
phases: 



o) = K'b)'''exp«), 



(-^b 



(A26) 
(A27) 



The oscillatory amplitude A(fevoive), which provides the only 
contribution to the remnant condensate density in Eq. ( IA25t 
with a significant dependence on fevoive, is then determined, in 
accordance with Eq. JA21> . by: 



A(fevolve) = «^«0<b)'^' ,/(^=.™,.+A,-./2) ^ (^28) 

Here the total phase shift 



A<p = 7T/2 + </)*'b + <Plfi + 2 [vih) - 'P(to)] 



(A29) 



consists of the sum of the phase shifts associated with the 
two-body amplitudes in Eqs. ( IA27t and ( IA26> . and twice the 
phase shift of the atomic mean field accumulated during the 
first magnetic field pulse of the sequence. This can be written 

as 



2[(p{h)-ip{to)] =2jui(fi -to)/ti 

-21m[Ai(fi -to) + Di(ti 



to)]. (A30) 



We note that the amplitude A(fevoive) on the right hand side of 
Eq. JA28I I depends just on the density of the gas rather than 

^bo^ob) inde- 
pendent of 

The angular frequency of the fringes is determined by 



\/n + 2[>f{t2)->p{ti)] litn-h). 



(A31) 



""V are the total number where E\ 



'evolve 



1^ is the binding energy of the highest excited di- 
atomic vibrational molecular bound state, and ^(fi) - <f(t\) is 
the phase shift of the atomic mean field accumulated during 
the evolution period. The assumption of weak interactions be- 
tween the atoms during the evolution period implies that the 
frequency shift due to the atomic mean field can be estimated 
in terms of the mean field chemical potential, i.e. by 

2 [ifiito) - ipih)] /(f2 - fl) ~ 2jU2/ri = 2jUevolve/?!. (A32) 



This estimate presupposes that the contributions of Eqs. ( IA15t 
and ( IA16> to the phase shift have decayed during the evolution 
period. This stationary limit, however, is never quite reached 
in the experiments J8|] . 
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6. Visibility of tlie Ramsey fringes in tlie remnant condensate 

Neglecting the dependence of the diffusion term Z)(fevoive) 
on the evolution time fevoive in Eq. ( IA25> . the remnant conden- 
sate density is an oscillatory function of fevoive with maxima 
and minima at evolution times f™'"; and r™", , which satisfy 

p vnlvp evolve' J 

the conditions of Eqs. ( I16> and \Y1\ . respectively. In accor- 
dance with Eq. jA25> , the visibility of the Ramsey fringes in 
the remnant condensate is then given by 



amsey 



l-exp(-2Re|A(e;..)-A(C,„)|) 

)|)' 



1 -t- exp (-2Re [A( 
which leads directly to Eq. ( I27> . 



*max \ A( *min 



(A33) 
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APPENDIX B: TWO-BODY TRANSITION AMPLITUDES 

In this appendix we give an analytic treatment of the two- 
body transition amplitudes for molecular association and dis- 
sociation as well as for elastic scattering at constant mag- 
netic field strength. These amplitudes can be used to provide 
estimates of the fringe visibility and the coupling function 
h(t, t) of the dynamic equation j28t for the atomic mean field 
[cf. Eqs. ( I27> and (IA13> . respectively]. The estimates apply to 
an idealised pulse sequence consisting of periods of constant 
magnetic field strength with sudden changes in between. 



1. Amplitude for elastic collisions 

We first consider elastic scattering because the derivation 
is more concise than for molecular association or dissocia- 
tion, while the basics issues are the same. In accordance with 
Eq. J21l i. the elastic transition amplitude is given by 



«1 dThi{T). (Bl) 

This result for a constant magnetic field strength B, (/ - 
1,2,3) of an idealised pulse sequence is well approximated 
in terms of a coupling function /i,(t) [cf. Eq. ( I20H in the con- 
tinuum limit. The derivation of Eq. JBll i takes advantage of 
the general formula 

+ \ f'dtUo(ti-t)V(Bi)Ui(t-ti^i), (B2) 



where Uoit, - f,-i) is the free time evolution operator. The 
latter can be expressed in terms of its spectral decomposition, 
i.e. 



t/o(f/-f/-i) = dp\p}e 



<pl- (B3) 



Here {r|p) = exp(/p ■ r/K)/(2nfi)^^^ is a plane wave associated 
with the relative momentum p and with the spatial relative 
coordinates r of an atom pair. 

Since the magnetic field strength is assumed to be constant 
in time, the time evolution operator t/,(f - f,_i) on the right 
hand side of Eq. (IB2t is given by the spectral representation 
of Eq. ( IA5> , which determines the coupling function to be: 



hi(T) = {inny 



I 



dp <0|y(fi,)l<\fiO)e 



{Q\V(Bd\MBd)e 



-iEb{Bi)T/n, 



b(B,)|0> 



(B4) 



Here we have substituted t - f,_i = t, which corresponds to 
the integration variable on the right hand side of Eq. (IB It . We 
use the Lippmann-Schwinger equation j44ll 

\4'p\Bd} =IP) + Go(p^/m + iO)V(Bd\(f>'^^\Bd) 

=|p) + Go(p^/m + iO)T(p^/m + /0)|p) (B5) 

to relate the continuum energy states of the relative motion of 
an atom pair (p^p\Bi) to the energy dependent T matrix. Here 



Go(z) ^{z + n^V^/m 



dp |p>- 



(Pi (B6) 



is the interaction free Green's function and "z - p~/m + /O" 
indicates that the energy argument approaches the kinetic en- 
ergy p^/m from the upper half of the complex plane. The 
choice of this particular set of continuum energy states is 
consistent with the asymptotic behaviour of their wave func- 
tions, at large distances r between the atoms, as described by 
Eq. dX^ . 

The plane wave matrix elements of the energy dependent T 
matrix can be expanded in terms of the parameter pa{Bi)/fi, 
which to first order accuracy leads to the universal self consis- 
tent approximation 



(27rS)^(Pout|rOr/m + /0)|pi„> 



47rfp-a{Bi)/m 
1 + ipa{Bi)lti' 



(B7) 



We note that the incoming and outgoing momenta pi,, and Pout, 
respectively, do not affect the approximate result even when 
their associated energies p^Jrn and pl^^Jm do not match the 
energy argument p~/m of the T matrix. Equation JB7l i is usu- 
ally referred to as the contact potential approximation. It re- 
flects, however, the universal low energy properties common 
to all T matrices, independent of the particular shape of the 
short range potential. 

Inserting Eqs. iB5\ and ( IB7> into Eq. ( IB4> . the right hand 
side of Eq. (IBU can be represented in the form: 

T2b(4>0 ^ 0O) =1 + Ai(ti - f,_i) + Di(ti - f;_i) 

-ifii(ti-ti^,)/n. (B8) 

The last contribution on the right hand side of Eq. jB8t stems 
from the plane wave part of the scattering wave function in 
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Eq. iB5\ and recovers the chemical potential ju,- of Eq. JA14> . 
via the general relation 

AnhMBd/m = (2nn)\0\V(Bi)\<t>^^^\Bi)), (B9) 

when the volume ^ is identified with the inverse condensate 
density (cf. SectionHli. 



The remaining part of the momentum integral over con- 
tinuum states on the right hand side of Eq. JB4> involves 
complex Gaussian integrals and determines the diffusion term 
D,(f, - f,_i). Expressed in terms of the complex error function 
of Eq. ( IA17> . the diffusion term is given by 



f>^ Jo J P /m ' ' 

iAntiaiBd j^'-'- ^ ^ . ^fj^^^^g^ (BIO) 



rrfV 



in the contact potential approximation. Here we have used 
Eq. iB5l as well as the universal relation 



£b(5,) = -h^/ [ma\Bi) 



(Bll) 



between the binding energy and the scattering length in the 
contact potential approximation (cf. Fig. The remaining 
time integral on the right hand side of Eq. ( IBlOt is readily 
performed using 



(B12) 



J"^ dy w (^i ^/iyj - i 1-2 y — - w {i yfix^ 



which leads directly to Eq. ( IA16> . 

The bound state contribution of Eq. ( IB8> is found from 
Eq. (|B4|i to be 



-;£b(B,)(f,-/,_,)/fi 



(B13) 



where we have used the momentum space Schrodinger equa- 
tion 

<p|y(B,)|<^b(B,)> = [£b(B,) - /^'H <p|0b(B,)) (B14) 

of the bound state wave function. Inserting the universal for- 
mula 



(r) = 



1 



^lna{B^ r 



(B15) 



for the bound state wave function in the contact potential ap- 
proximation into Eq. JB13> yields Eq. ( IA15t . 

We note that in applications to time independent low en- 
ergy two-body collision phenomena the contact potential ap- 
proximation is valid only in a small region of magnetic field 
strengths in the close vicinity of the zero energy resonance, as 
indicated in Fig.|3l In this region the scattering length exceeds 
by far all the other length scales set by the inter-atomic inter- 
actions. A second restriction of the contact potential approx- 
imation, specific to the considerations in this paper, is related 
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FIG. 1 1 : Comparison between the contact potential approximation 
and exact numerical calculations of the two-body time evolution op- 
erator flolfl^EcIl in applications to the first order iterative approxi- 
mation of the atomic mean field [cf. Eq. jAlOH . The numerical 
calculations account for the length scales set by the scattering length 
and by the van der Waals interaction (cf. Fig.|3}. The initial conden- 
sate density of «c('o) = 3.16 x 10'^ cm"' is obtained from the low 
density experiments of Ref. 01, while the stationary magnetic field 
strength of 6i = 155.5 G corresponds to the typical minimum of the 
first pulse of the sequence in Fig.Q 



to the fact that dynamical phenomena involve a whole range of 
inter-particle collision energies E. The approximation is ap- 
plicable only at low collision energies, which depend on the 
magnitude of the scattering length through pm^^a{Bi)IJi «: 1. 
Here p^^^ - V»4£maxJ denotes an estimate of the maximum 
momentum scale set by the time evolution of the observed 
quantity. Figure ^2 compares the contact potenti al ap proxi- 
mation with exact numerical calculations llOt [l9l 12011 of the 
two-body time evolution in applications to the first order it- 
erative approximation of the atomic mean field ^^(f) (cf. Ap- 
pendix|XJ. The comparison confirms that the approximations 
leading to the contact potential approach are well satisfied for 
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the elastic transition matrix element of Eq. iBl\ at the typical 
magnetic field strength of Bi = 155.5 G during the first pulse 
of the sequence (cf. Fig.Q. 

2. Molecular production and dissociation in magnetic field 
pulses 

In this subsection we estimate the probabilities and 
P^Q, in the contact potential approximation, for idealised 
magnetic field pulses, which just involve sudden changes of 
the magnetic field strength. The derivations provide analytic 
estimates that reveal the sensitivity of the fringe visibilities 



[cf. Section |lll to the shape of the magnetic field pulses. 
We shall consider just the amplitude T^^{4>1™^^'' <- if>o) of 
the first magnetic field pulse, whose time evolution operator 
Ui{ti - to) = C/#i(fi,fo) is translationally invariant in time, 
under the assumption of an idealised sequence. The second 
amplitude T^{(p() «— 0^*°'^'^) can then be obtained from con- 
siderations of time reversal. 

As the determination of r|g((/i^™''''= «- 0o) closely follows 
the considerations in the preceding section, we shall just sum- 
marise briefly the main elements of the derivation. Taking 
advantage of the decomposition of the time evolution oper- 
ator U\{t\ - to) in Eq. JB2> . the transition amplitude in the 
continuum limit is given by: 



r#l / jL evolve , f \ 
2BWh <Po) 



(Inti)^ 



h(B2)m + 



in 



dT <0b(S2)|t/o(fl -to- T)V(Bl)Ul(T)\0} 



(B16) 



In analogy to the derivation of the amplitude for elastic col- 
lisions, we replace the time evolution operators U\{t) and 
Uo{t\ - to - t) by their spectral decompositions in Eqs. (IA5> 
and iB3\ . respectively. The energy dependent T matrix el- 
ements that result from the continuum part of the two-body 
spectrum can be approximated in complete analogy to the pre- 
ceding section. The momentum space bound state wave func- 
tion can be determined from the Schrodinger equation ( IB14> 
to be 



(<Pb(B)\p} ^ 



{UB)\V(B)\p} 
EbiB) - pVm 



EbiB) 



(B17) 



n{m\E^{B)\)ynE^{B)-pym] 
in the contact potential approximation. The resulting momen- 



tum and time integrals are similar to those in the preceding 
section. 

In the contact potential approximation, the transition ampli- 
tude of Eq. (IB16> can be expressed in terms of; the scattering 
length fl(Bi) of the first magnetic field pulse, the ratio of scat- 
tering lengths 



r = a(B2)la(B,) 



(B18) 



associated with the first magnetic field pulse and the evo- 
lution period, respectively, in addition to the phase shifts 
6i = |£b(Bi)l(fi - to)in and 62 = |£b(B2)l(ri - to)lti. We note 
that just three of these parameters can be independently var- 
ied. The final formula for the transition amplitude then reads 
[cf. Eq. JA17> and Ref. ll45ll for a discussion of the properties 
of the function w(z)]: 



7^#l / jevolve 



00) 



87ra3(fi,)y 



y - \ +w{i -\IJ62j + H 2 (' ~ ^ (' V^)] 

- 162 I dx w (i -\/i62{l - X)] - w ii -\/id2X 

Jo [ yi7T62x 



(B19) 



The application of the analytic formula of Eq. ( IB19> to the 

molecular production efficiency P^^^ - |7'2B'-^b™'™ *~ '/'o)|" 
in Fig. illustrates that the contact potential approximation 
follows the qualitative trends of the exact numerical result. 
Equation ( IB19> clearly reveals that the efficiency of molecu- 
lar production depends on the ratio fl(Bi)/['y] '^^ of the length 
scales a{Bi) set by the low energy inter-atomic interactions 
and the mean inter-atomic distance, respectively, on one hand 



(cf. Sectionlnlfor a discussion of the physical meaning of the 
volume in applications to dilute gases), and on the ratio 
y = a(B2)la{Bi) on the other hand. The scattering length 
a{B2) thereby provides an estimate of twice the mean internu- 
clear distance associated with the bound state wave 

function 0^^°'^'^ in the evolution period. An optimal molecular 
production efficiency thus requires a proper balance between 
these length scales, which can differ by orders of magnitude 
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• exact 

— contact potential 



^ 9 




under the experimental conditions of Refs. 



500 



FIG. 12: 



The 



p#i 



molecular production efficiency 

as a function of the hold time ti - tg of the 
first magnetic field pulse for a minimum magnetic field strength 
of 5min - 155.5 G (cf. Fig. The magnetic field strength in 
the evolution period is Seyoive = 156.5 G. The figure shows exact 
numerical calculations (filled circles) fiol l20ll and the contact 
potential approximation obtained from Eq. <B19t . We note that the 
product 'V/^'i, is independent of the volume 'V. 



We also note that Eq. ( IB19> is less accurate than Eq. (IB8> 
for the elastic amplitude. The reason for this is mainly related 
to the different energy scales that enter the different ampli- 
tudes. In the case of molecular production the maximum en- 
ergy scale is set by the binding energy \E\,(B2)\ = Pmax/"^ ™ 
the evolution period. The requirement pmiLxa{B[)/fi 1, as- 
sociated with the contact potential approximation, is violated 
in the applications. Studies beyond the scope of this paper 
show that this violation mainly affects the evolution on short 
time scales on the order of fi - fo ^ fi/\Eb{B2)\- Despite this 
deficit, however, the contact potential approximation provides 
useful qualitative insight into the dependence of the molec- 
ular production efficiency on the very different length scales 
set by the microscopic binary collision physics and the mean 
inter-atomic distance of the gas, respectively. 
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